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Black holes with masses ~ 1 Mo cannot be produced via stellar evolution. A popular scenario 
of their formation involves transmutation of neutron stars — by accumulation of dark matter trig- 
gering gravitational collapse in the star centers. We show that this scenario can be realized in the 
models of bosonic dark matter despite the apparently contradicting requirements on the interactions 
of dark matter particles: on the one hand, they should couple to neutrons strongly enough to be 
captured inside the neutron stars, on the other, their loop—induced self-interactions impede col- 
lapse. Observing that these conflicting conditions are imposed at different scales, we demonstrate 
that models with efficient accumulation of dark matter can be deformed at large fields to make 
unavoidable its subsequent collapse into a black hole. Workable examples include weakly coupled 


models with bended infinite valleys. 


I. INTRODUCTION 


The existing gravitational wave detectors are sensitive 
to compact objects — black holes (BHs) and neutron 
stars (NSs) — of masses ranging from tens of solar masses 
down to Mo, see [1, 2]. Thus far, few tens of coalesc- 
ing BHs and several BH—NS systems were detected [3-6], 
while hundreds of observations are further expected in 
the near future. This will be sufficient to map out the 
mass function of the stellar—size BHs. 

Standard stellar evolution predicts no BHs lighter than 
the heaviest neutron stars, as the Fermi pressure stabi- 
lizes the neutron cores at masses below 2.5Mo, cf. [7-9]. 
The BH mass function is hence expected to have a gap 
below this value, and any detection in that region would 
automatically imply existence of an exotic mechanism for 
black hole formation. We will refer to such light BHs as 
“solar mass black holes.” 

A widely discussed possibility is production of solar 
mass BHs in the early Universe from small-scale density 
perturbations. It has a natural extension: these objects 
may also play the role of dark matter (DM). The lat- 
ter hypothesis, however, is strongly constrained by lens- 
ing [10], CMB measurements [11, 12], and even the grav- 
itational wave observations themselves [13]. 

Alternatively, the solar mass BHs may appear in the 
present Universe due to neutron star collapses caused by 
seed BHs [14]. The latter BHs may in turn be of primor- 
dial origin, and be small enough to evade the lensing and 
CMB constraints. This option, however, still relies on 
cosmology to produce a sufficient seed BH population. 

A more intriguing possibility is the formation of seed 
BHs through the gravitational collapse of dark matter ac- 
cumulated by old neutron stars [15-21]. A sufficient ac- 


cumulation is possible only if the DM is asymmetric [22] 
and does not annihilate in the stellar core. In the simplest 
case one assumes that the DM particles carry a conserved 
global charge ensuring their stability, and their antipar- 
ticles disappear during the early stages of cosmological 
evolution. 


In the present paper we focus on the last scenario. Our 
goal is to identify DM models where such catalyzed trans- 
mutation of the neutron stars into the solar mass BHs can 
take place. The mechanism we consider involves four 
stages: (i) capture of DM by the neutron star; (ii) its 
thermalization; (iii) concentration in the star center and 
gravitational collapse into a seed BH; (iv) accretion of 
the star material onto the BH and formation of the solar 
mass BH. Given the DM abundance, the stage (i) sets the 
maximum mass of DM that can be accumulated by the 
neutron star. On the other hand, the stage (iii) requires 
some minimum amount of DM for the successful grav- 
itational collapse. The viability of the overall scenario 
depends on the possibility to satisfy both requirements 
in the same DM model. 

Stages (i), (ii), and (iv), ie. the beginning of the pro- 
cess and its very end, have been extensively discussed 
in the literature, see the review [23] and [24, 25]. We 
will provide a short summary of these results in the next 
section. The main part of this paper considers model- 
dependent evolution and gravitational collapse of dense 
dark matter cloud at stage (iii). 

The minimum number of particles required for col- 
lapse is different in the cases of bosonic and fermionic 
dark matter, and crucially depends on the DM self- 
interactions. In the ensemble of free DM fermions, the 
Fermi pressure balances self-gravity and halts the gravi- 
tational collapse unless the total particle number exceeds 


the (analog of the) Chandrasekhar limit [26], 
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where m is the DM mass. This number is prohibitively 
large unless m is in the 100 TeV range or above [19]. 
While not necessarily impossible, such scenarios with 
super-heavy asymmetric dark matter are outside of the 
mainstream cosmological models, and we do not consider 
them here. 

If the DM is bosonic and non-interacting, the grav- 
itational attraction of its particles is balanced by the 
kinetic (“quantum”) pressure guaranteed by the uncer- 
tainty principle. Gravitational collapse happens if the 
particle number exceeds the critical value 
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which is parametrically smaller than in the case of 
fermions. Bosonic DM is therefore a more promising can- 
didate for the creation of seed BHs. In the rest of this 
paper we will focus on bosons. 

While looking advantageous in the free case, Eq. (2), 
however, can be ruined by very weak DM self- 
interactions which are generically present in our scenario. 
Indeed, DM capture and thermalization at stages (i) and 
(ii) are possible only for sufficiently strong DM-nucleon 
interactions which in turn induce DM self-couplings via 
loops [18]. In the simple one-field models, the induced 
terms in the DM potential are quartic in fields and should 
be positive, i.e. repulsive, for vacuum stability. In ef- 
fect, they generically oppose the gravitational collapse 
and increase the required critical DM multiplicity to un- 
acceptably large values. We have therefore conflicting 
conditions on the DM model. 

In this paper we show that the conflict can be resolved. 
Namely, the number of DM particles required for collapse 
can be almost as small as in Eq. (2), and all of these 
particles can be accumulated inside the neutron star by 
scattering off nucleons. We start by clarifying and quan- 
tifying the requirements on the DM models. Then we 
propose a generic mechanism to satisfy them simulta- 
neously. We find that the smallest DM multiplicity for 
collapse is achieved in the models where: (i) the DM po- 
tential includes a long valley extending to large, in some 
cases Planckian fields, and the cutoff of this potential 
is high enough; (ii) both attractive and repulsive self- 
interactions are suppressed along this valley; (iii) loop 
quantum corrections from interactions with the visible 
sector do not break the condition (ii). 

Our mechanism is based on the fact that the conflict- 
ing requirements are imposed at different scales. The DM 
capture depends on the physics in the vicinity of vacuum, 
while self-interactions obstructing the gravitational col- 
lapse should vanish at strong fields. Thus, one can make 
the potential valley bended in the field space. It may start 
going along the dark matter field ¢, then take a turn at 


some |¢| Z A and continue in the direction of another 
field x. In this case the interaction of ¢-particles with 
nucleons, as we will argue, does not generate an effective 
potential for x, and the latter may grow and collapse into 
a BH almost as in the free bosonic case. 

The rest of this paper is organized as follows. We re- 
view accumulation and thermalization of DM inside the 
neutron stars in Sec. II. General requirements on the DM 
models with the smallest critical multiplicities for grav- 
itational collapse are derived in Sec. III. In Sec IV we 
propose a mechanism to satisfy these requirements. We 
conclude in Sec. V. 


II. DM INSIDE THE NEUTRON STAR 


In this section we summarize the existing results on 
DM capture and thermalization in the neutron star. Two 
points will be essential for us: (i) the total amount of 
captured DM and its dependence on the strength of DM- 
neutron interactions; (ii) the fact that with time this DM 
forms a Bose-Einstein condensate described by classical 
fields. 

For concreteness we will assume that the DM particles 
are globally charged scalars with mass m belonging to 
the typical WIMP range from GeV to a few TeV. 


A. DM capture 


Two generic mechanisms trap DM inside the neutron 
star: accumulation during the star lifetime and gravita- 
tional capture at the star formation. These mechanisms 
are cumulative. Potentially, they may provide compara- 
ble amounts of trapped DM. 

In the next sections we consider accumulation of DM 
particles from the ambient galactic halo [27, 28]. Far 
away from the neutron star, the distribution of DM ve- 
locities is nearly Maxwellian with small dispersion 07, 
e.g. U ~ 107? in the Milky Way. But the particles ac- 
quire semi-relativistic speeds ~ 0.5 as they fall into the 
neutron star. With masses in the WIMP range, they 
lose energies of order Mp — the neutron’s rest mass — 
in collisions with neutrons. This is much larger than the 
asymptotic energies of the particles; thus, most of them 
bind gravitationally to the neutron star after the first 
collision. Besides, the momentum transfer in their col- 
lisions marginally exceeds the Fermi momentum of neu- 
trons, so that the degeneracy of the latter does not play 
a crucial role. Neglecting the degeneracy is a crude ap- 
proximation that overestimates the capture rate [29-31], 
but in view of comparable astrophysical uncertainties we 
will use it for simplicity. Finally, we will ignore general 
relativity effects which enhance the capture rate by an 
order 1 number, see Ref. [32]. 

Under these assumptions, the capture rate, i.e. the 
number of trapped DM particles per unit time, takes the 


form [27, 28], 


dN 
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Here ppm is ambient DM density, M, and R, represent 
the mass and the radius of the neutron star, and we in- 
troduced the probability f < 1 for the DM particle to 
scatter during one pass through the star. The value of f 
is proportional to the DM-—neutron cross section o, 

f = O/C crit if O < Ocrit 5 (4) 
and f = 1 otherwise. Here the proportionality coefficient 
(“critical” cross section [27]) depends on the neutron star 
parameters, 


Gorit ~ R2m,/M, ~ 107% cm? , 


where we used M, ~ 1.5 Mo and R, ~ 10 km in the es- 
timate. Notably, Gerit is comparable to the upper limits 
on the DM-nucleon cross section coming from the direct 
detection experiments [33-37]. 

Since dN/dt x mt, the accumulation rate of dark 
mass does not directly depend on m. Integrating Eq. (3), 
we obtain the total DM mass inside the 10 Gyr old neu- 
tron star, 


MW: Miot ~ 7 x 10° GeVf ~ 7x 10°? Mof, (5) 
dwarf: Miot ~ 8 x 1046 GeV f ~ 7 x 1071 Mof . (6) 


The above two options differ by the parameters of 
the ambient DM distribution: t =~ 220 km/s, 
pom ~ 0.3 GeV/cm? in the Milky Way (MW) and 
ü~ 7 km/s, ppm ~ 100 GeV/cm? in the densest dwarf 
galaxies [38, 39]. Besides, taking f = 1 turns the above 
estimates into model-independent upper limits on the to- 
tal DM mass that can be captured by the neutron star 
in the respective environments. 

Another possibility is a gravitational trapping of the 
DM during formation of the neutron star [40, 41]. The 
latter objects are created in the supernova collapses of or- 
dinary stars with masses above ~ 9Mo, which, in turn, 
are born in giant molecular clouds. Originally, as the 
baryonic gas contracts adiabatically into a proca star, 
a low-—velocity fraction of the ambient DM gets trapped 
by its gravitational well. Eventually, this DM ends up in 
the center of a heavy progenitor star; estimates of [40, 41] 
show that the captured DM mass is close to Eqs. (5), (6) 
within an order of magnitude. Later, a fraction of the 
star DM is inherited by the neutron star in the course of 
supernova collapse. That last process has never been con- 
sidered in the literature, but our rough estimate suggests 
that the respective suppression of Mi 4 lies between 1 
and 107°. 

In the rest of this paper we will use Eqs. (5), (6) 
and ignore the second, “gravitational capture” mecha- 
nism despite the fact that it seems less sensitive to non- 
gravitational DM interactions. Indeed, that additional 
mechanism cannot dramatically change the amount of 


captured DM even in the best case. Besides, it de- 
pends on the multi-stage neutron star formation which 
is subject to large astrophysical uncertainties. Finally, it 
does in fact rely on the DM—DM and DM-neutron in- 
teractions, albeit in a subtle and indirect way: the non- 
gravitational couplings are needed to thermalize the DM 
particles inside the progenitor star and mix them in the 
phase space, or none would lose enough energy to get into 
the neutron star. 


B. Thermalization and condensation 


Once gravitationally bound, the DM particles settle 
on star—crossing orbits and continue to lose energies in 
repeating collisions with neutrons until a thermal equilib- 
rium is reached. Their kinetic energies reduce to the neu- 
tron star temperature T ~ 10° K and their orbits shrink 
to the “thermal” radius, 


T 1/2 1 1/2 
Tih © ( : ) ~ 20cm (=s) » 7) 
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where we substituted the density pe ~ 1015 g/cm? of the 
neutron star core. The particles steadily arrive into the 
central DM cloud during the entire lifetime of the neutron 
star with the rate (3). 

The characteristic time of this thermalization process 
was estimated in Refs. [14, 42]. For almost any model 
within our scenario, it is much shorter than the age of the 
Universe [14]. Indeed, we had already mentioned that 
BH formation requires a large number of DM particles 
inside the neutron star [18]. To capture all of them, one 
usually assumes the largest possible DM-neutron cross 
section 0 ~ Cerit, and even that may be insufficient. With 
these interactions, the DM particles equilibrate quickly, 
since every star crossing in the beginning of the process 
leads to scattering. On the other hand, we will be able to 
consider very small DM-neutron cross sections once the 
new mechanism for reducing the required multiplicity to 
Eq. (2) is invoked. In that case a necessity to thermalize 
the DM imposes the strongest constraint on its interac- 
tions with neutrons. 

As the DM particles continue to accumulate, the mul- 
tiplicity N of the central cloud grows at a fixed ra- 
dius (7). Eventually, the mean distance between the par- 
ticles rm N13 drops below the size (mT)~!/? of their 
wave functions. This happens at 


BEC: N 2 0.2 T’ M pz”? ~ 10%, (8) 
i.e. significantly before the maximal number of parti- 
cles (5), (6) is reached. At this point the particle wave 
functions start to overlap and Bose-Einstein condensate 
forms in the thermal cloud [43-46]. Since then, most 
of the thermalized DM particles occupy the lowest level 
in the combined DM and neutron star gravitational po- 
tential, with the thermal energy being carried by a few 
remaining particles. 


Thanks to large occupation numbers (8), the Bose- 
Einstein condensate at the lowest level can be described 
by a classical DM field ¢(x). From the very beginning 
it forms a non-rotating [47] self-bound soliton which is 
almost detached from the neutron star surroundings. In- 
deed, even at multiplicity (8) the gravitational field of 
this object can be estimated to exceed the neutron star’s, 
and with time the soliton mass grows. We will call 
this soliton a Bose star [48, 49] or a Q-ball [50-52] if 
it is mostly bound by self-gravity or by attractive self- 
interactions, respectively. In either case, the properties 
of the soliton can be determined by solving the stationary 
classical field equations, where a single free parameter is 
the total number of accumulated DM particles N. 

With the soliton formation the final, model—dependent 
and nonlinear, stage of DM evolution begins. As N in- 
creases from Eq. (8) to the maximal values in Eqs. (5), 
(6), the soliton becomes heavier. It collapses gravitation- 
ally if the condition for collapse — the hoop conjecture — 
gets satisfied prior to accumulating the maximal amount 
of DM. If, to the contrary, the soliton size exceeds its 
Schwarzschild radius even for the largest N in Eq. (6), 
the black hole does not form. 

Let us finish Section by commenting on DM thermal- 
ization in white dwarfs — second to best compact ob- 
jects accumulating DM. Their cores are significantly more 
dilute, pe < 108 g/cm, and have higher typical tem- 
peratures T > 10°K as compared to the neutron stars, 
see [53]. Thus, Bose-Einstein condensation of dark mat- 
ter particles in their centers would require larger DM 
multiplicity N > 104°, see Eq. (8). But in fact, the 
DM inside the thermal radius becomes self-gravitating 
before that, at N > 104° (m/100 GeV)~*/?. It collapses 
gravithermally and forms a compact object where the 
subsequent DM cooling and (possibly) Bose-Einstein 
condensation occur. This process deserves a separate 
study, which is outside of the scope of this paper. 


III. SELECTION RULES 
A. Optimizing the DM model 


As we stressed in the Introduction, interactions gener- 
ically detain the gravitational collapse until the number 
of DM particles becomes parametrically larger than in 
Eq. (2), and this is often too large to be accumulated by 
the neutron star during the Universe’s lifetime. To gain 
a more quantitative understanding of the numbers, we 
start with the simplest scalar DM described by a single 
complex field, 


£L=|V, 9)? —V(\¢)). (9) 


We will try to optimize its scalar potential V in a way 
that minimizes the dark matter multiplicity required for 
collapse. We consider minimal coupling to gravity and 
ignore interactions with the visible matter — those will 


(a) attraction (b) repulsion 
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FIG. 1. Possible forms of the scalar potential. 


be added later. Importantly, we also assume that the 
typical scale(s) A in the potential V are essentially sub- 
Planckian, 


This separates our model from the effects of quantum 
gravity. 

Suppose a compact, stationary, and stable solitonic 
configuration of the scalar field — a Q-ball or a Bose 
star — is formed inside the neutron star. Let us com- 
pute its mass M and radius R as functions of the global 
charge N. By conservation law, the latter quantity 
counts the number of dark matter particles appended 
to the Bose-Einstein condensate. We recall, first, that a 
generic stationary solution in the model (9) has the form, 


= (x) e, (11) 


Here we introduced the real-valued soliton profile y(a) 
and the energy w of particles inside it. Second, the 
gravitational field of the soliton is expected to be small, 
Juv © Nv, except for the critical case when gravitational 
collapse is about to happen. In this case the flat-space 
Noether charge can be used: 


Nw if ëa (b* ob — b00¢*) ~ wpe R? , (12) 


where the last equality is a crude estimate in terms of 
the soliton size R and the field yo ~ (0) in its center. 
Since every particle inside the soliton has energy w, the 
total mass of this object is of order M ~ wN. 

Recall, however, that the soliton parameters yo and 
R are related by the field equations which involve the 
scalar potential V (po). We therefore consider two op- 
tions. First, one can assume that the potential grows 
almost quadratically at weak fields, V ~ m?|¢|?, and 
then flattens out beyond some scale Aatt like V œ |d|* 
with a < 2, see Fig. l(a). This case corresponds to parti- 
cle attraction inside the solitonic core, as the energy per 
unit charge is smaller than the particle mass. At strong 
fields we approximate, 


Vw m?o Ase? at o |b] Z Aas (13) 


where 0 < a < 2. It is precisely the scalar self—interaction 
that holds the soliton — the Coleman’s Q-ball [51] — 
together, since gravity is weaker: Aatt < M, according 
to Eq. (10). In the field equation, the self—attraction is 
balanced by the kinetic pressure: 0?¢ ~ V’ ~ w?¢, or 


po/R? ~ mpg" Man” ~ w eo, (14) 
see Appendix A for details. This gives w ~ R71, 


N ~ (RA)? (Rm) C9, and M~ N/R, (15) 


where we parametrize the soliton configurations with 
their sizes R. 

The soliton collapses gravitationally when its compact- 
ness R,/R = 2M/(M;,R) becomes of order one, i.e. at 
masses above critical, 


Mer © RM},/2 ~ wNer - (16) 


Substituting this relation into Eq. (15), we obtain the 
number of DM particles needed for collapse, 

M2, M l 2-a 
= (=) for 0<a<2. (17) 
At the critical point, the field inside the soliton is Planck- 
ian: Yo,er ~ Mp, cf. Eqs. (14), (15), and (17). 

Notably, the critical multiplicity (17) is minimal at 
a = 2 when the “free bosonic” expression (2) is recov- 
ered. The other values of a < 2 are less advantageous, as 
Aatt is parametrically below the Planck scale. One ob- 
tains almost the “fermionic” multiplicity (1) in the case 
V œ |d| and even larger Ner for the -independent po- 
tential. Thus, contrary to naive expectations particle 
attraction obstructs collapse, the reason being that the 
energy per particle w becomes much smaller than m. 

One can push the above “attractive” option to the ex- 
treme assuming that the scalar potential V (¢) decreases 
at strong fields, e.g. V x —|¢|*. However, that would 
destabilize the soliton making its field evolve towards 
lower V, cf. [54, 55]. With no new positive terms in the 
potential to stop the process, the region with V (¢) < 0 
would be eventually reached [56]; then the soliton turns 
into an expanding and Universe-destroying bubble of 
true vacuum [57]. Even if the bubble can be somehow 
forced to collapse gravitationally, the region in its cen- 
ter still breaks the positivity conditions, so that a naked 
singularity may appear instead of a black hole. If, al- 
ternatively, the potential starts growing at larger fields, 
again, the field stops rolling at that point, thus bringing 
us back to the two options in Fig. 1. 

Equation (17) hints at the possibility that the second 
option of repulsive self—interactions with a > 2 may be 
more interesting, see Fig. 1(b). In this case the only at- 
tractive force is gravity. The respective soliton is called 
a Bose star [48, 49], since it is bound by gravitational at- 
traction compensating interaction pressure of Bose parti- 
cles inside it. Note that the repulsive energy gives a sub- 
dominant contribution to the mass of the subcritical ob- 
ject because its opponent — the gravitational energy — 


Ner m 


remains small until the rim of collapse. We therefore keep 
two terms in the potential, 
V & mg)? + m?|6|* Ar (18) 


rep ? 


where now a > 2, the field ¢ is arbitrary, and Arep sat- 
isfies Eq. (10). Performing the estimates similar to the 
ones before [58], we find out that the Bose star collapses 
gravitationally only at a > 8/3, see Appendix A for de- 
tails. In this case the soliton field indeed gets stuck in the 
region of subdominant self-interactions yo S Arep until 
the collapse, at which point Yo, cer ~ Arep and the soliton 
charge equals 


pi for a > 8/3. (19) 


Notably, this critical multiplicity is independent of a, see 
also Refs. [59] and [58]. At Apep ~ m the expression (19) 
reproduces the fermionic result (1). The way to decrease 
the multiplicity is to increase Arep suppressing the self- 
repulsion. At Arep ~ Mp; that force is as weak as gravity, 
the field inside the critical Bose star is Planckian, and we 
obtain the “free bosonic” formula (2), again. 

The remaining region 2 < a < 8/3 is worst of them 
all, since in that case the Bose-Einstein condensate does 
not clump under self—gravity. Indeed, an estimate of Ap- 
pendix A shows that the self-repulsion |¢|* with a in 
this range is stronger than gravity at large distances. It 
makes the DM condensate spread over the entire volume 
available inside the neutron star gravitational field. 

We conclude that the critical particle number Ner 
is smallest in the restricted class of models with long, 
Planckian-size valleys |¢| £ Mpı and almost quadratic 
potentials V œ |d|? at their bottoms. Any interaction 
impedes collapse and sharply increases the critical mul- 
tiplicity. In particular, self-repulsion becomes strong, al- 
most equivalent to the Fermi pressure at large fields — 
hence the “fermionic” result (19). Self-attraction does 
not help either: it provides negative binding energy and 
lowers the soliton mass which is bad for collapse. The 
respective critical particle number is also parametrically 
larger, see Eq. (17). 

In the optimal model with exactly quadratic poten- 
tial, the kinetic pressure inside the Bose star is com- 
pensated by the gravitational attraction [48, 60, 61]. 
The respective object becomes gravitationally unstable 
at Yo, er ~ Mpi- It has the critical multiplicity [60] 


Ner ~ 0.653 MZ,/m? fora=2. (20) 


This result agrees with the estimate in Eq. (2). Non- 
relativistic approximation remains valid during the most 
part of the Bose star growth and gets marginally broken 
at the rim of collapse. 


B. Obstruction by quantum corrections 


How far the DM model can deviate from the free 
bosonic theory? To find out, we require that the critical 


multiplicity for collapse Ner does not exceed the maximal 
amount (6) of captured DM: mMNer L Mawat, 

In the attractive case, this condition bounds from be- 
low the scale |¢| ~ Aatt at which the scalar potential 
flattens out in Fig. 1(a): 


2x 10-9 Gey \ C7 
BE i (21) 


see Eqs. (6), (17) and recall that 0 < œa < 2. Typi- 
cally, Aatt is very large. Indeed, since the critical Q- 
ball has yo ~ Mp, the model with flat potential should 
be trustable, i.e. renormalizable and weakly coupled, all 
the way up to the Planckian scale. The only manifestly 
renormalizable flat potential is V = const, it appears [62] 
e.g. in the celebrated Friedberg—Lee-Sirlin model [50]. In 
this case 


Natt ea My ( 


Aatt Z 5 x 1013 GeV (mf /100GeV)-/?,, a=0. (22) 
Soon we will see that large scales (21), (22) are problem- 
atic because multiloop corrections to the scalar potential 
become relevant at yo < Aatt. They generate interaction 
pressure and prevent collapse. 

The scale Agi, can be substantially lowered in a spe- 
cific class of renormalizable multifield models where the 
potentials at the bottoms of curved valleys have a 7 2. 
We will consider this option in Sec. IV. 

In the opposite, self-repulsive case the inequality 
mNer L MS“ also provides a large scale 


Arep Z 2 X 108 GeV (mf/100GeV)~*, (23) 


following from Eqs. (6) and (19). This condition 
strongly suppresses all repulsive self-interactions at fields 
lol S Arep, see Eq. (18). For example, in the \4|¢|4/4 
case, 


Ma = (2m/Arep)? S 7 x 10713 f2(m/100 GeV)*. (24) 


It is worth remarking that the renormalizability of the 
repulsive potential is not required, since the Bose star 
field yo S Arep remains parametrically below the natural 
cutoff of the theory (18) prior to collapse. 

Now, recall that our dark matter should interact with 
the visible sector, and strongly enough, in order to be 
captured by the neutron star. The respective couplings 
should be renormalizable and stay under control at strong 
fields — hence, their forms are constrained by the Stan- 
dard Model symmetries. Let us demonstrate that, gener- 
ically, loop corrections from these interactions break the 
desired properties of the dark matter potential: lift its 
flat parts and generate unacceptably large repulsive ver- 
tices, as has been first pointed out in [18]. 

Couple, e.g., the field ¢ to the Higgs doublet H (x) by 
deforming the potential of the latter to 
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FIG. 2. Scattering ọn —> ọn. 


cf. [63, 64]. Here v ~ 246 GeV and Ay ~% 0.13 are the 
standard Higgs parameters, and the new constant y reg- 
ulates its interactions with the dark sector. Notably, the 
model is weakly coupled at y < (4t\yq)!/? ~ 1. 

The physical Higgs field h(a) is defined as 
HtH =(v+h)?/2. It interacts with neutrons n 
via the effective Yukawa vertex Vann © 2m,hiin/(9v) 
which is known up to light quark contributions giving 
a factor 1 + 2 [65, 66]; m, œ~ GeV is the neutron 
mass. This means that the dark matter also scatters off 
neutrons. The respective diagram is shown in Fig. 2 and 
the cross section at nonrelativistic momenta equals! 

Qo 4 
pene i ae (26) 


~ 81ami,m? ’ 


where my = v(2\q)'/? ~ 125GeV. Importantly, the 
scattering probability f = o/oait < 1 should be large 
enough to capture the ¢—particles inside the neutron star, 
see the constraints (21), (23). Thus, the coupling y is 
large. 

On the other hand, the same coupling (25) generates 
DM pressure at strong fields. In fact, we have already 
tuned this potential to cancel DM self—interactions at the 
classical level. Namely, at large fixed yo = |¢| the Higgs 
field adjusts itself to minimize the first term in Eq. (25), 


v? yll? 
2° Qe 


1 
H'H = zh v)? = 


(27) 


As a consequence, the tree potential is quadratic at the 
bottom of this potential valley, V ~ m?|ø|?. Neverthe- 
less, the self-coupling reappears, again, once loop con- 
tributions from the visible matter are included. As an 
illustration, let us first ignore all fields except for the 
Higgs boson and ¢ itself. Then their one-loop effective 
potential [68] at large |ġ| « h along the valley (27) equals, 


AVi —loop x A4, eff |e|*/4 . (28) 


where Aq ef = y?L/(277) + O(y?/Az) is a correction to 
the |@|* self-coupling, L = In(|¢|/Ayen), and Aren ~ m is 
a renormalization scale. Note that a proper calculation of 


1 Hadronic form factors and strong nucleon interactions may sup- 
press this cross section by an additional factor of up to 1073, 
see [31, 67]. However, that would only sharpen the arguments of 
this section. 


the effective potential at strong fields includes renormal- 
ization group resummation of the leading logs [69, 70]. 
However, that usually introduces an order—one factor in 
Aa, eff Which does not affect our estimates. 

We see that Eq. (28) brings in the |¢|* repulsion even 
if it was absent before. Moreover, since A4, eg logarithmi- 
cally depends on the field, the new contribution cannot 
be canceled by any renormalizable counter-terms even if 
an arbitrarily precise fine-tuning is allowed. An obvi- 
ous way out is to make the one-loop repulsion small, so 
that it does not preclude black hole formation. Requir- 
ing A4, ef to satisfy Eq. (24) with o given by Eq. (26), we 
obtain the inequality 


y > 400 L/? > 400, (29) 


which cannot be satisfied in a weakly coupled model with 
y <1. Moreover, even this unacceptable lower limit can 
be achieved only if m = PeV. We conclude that either 
the coupling constant y is too small for accumulating 
the required amount of ¢—particles, or the interaction 
pressure caused by the same constant prevents the ¢-— 
condensate from collapsing. 

As an alternative, one may try to couple ¢ to fermions 
that generate negative terms in the effective potential. 
In the model (25) this amounts to recalling that every 
massive Standard Model field gives a potential to Higgs 
and therefore produces a vertex |@|* œ h* along the val- 
ley (27). The leading contribution is negative and comes 
from the top quark Yukawa coupling ye + 1. We obtain 
ON, of = —3y7y?L/(877AZ,). This definitely destabilizes 
the valley [71, 72] unless even larger positive |¢|* terms 
are introduced, which, however, returns us to the no-go 
estimates given above. 

Somewhat more elegantly, one may organize a (par- 
tial) cancellation between the fermionic and bosonic 
loops. But that would mean upgrading the Stan- 
dard Model to (N)MSSM. In that case the inequal- 
ity mN., < Mf imposes strong constraints on the 
supersymmetry—breaking operators that detune the can- 
cellation [18]. In Sec. IV we consider more economic and 
general possibility. 

To sum up, loop corrections from the visible sector are 
dangerous, generic, and cannot be avoided. Thus, we 
need a special mechanism to tame them. 


C. Requirements for the DM model 


Let us summarize the properties of dark matter model 
needed for black hole formation inside the neutron star. 

(a) The scalar potential of the model should include a 
long valley parametrized by the complex scalar ¢. The 
valley should extend to large fields: to |¢| ~ Mp in the 
case of attractive self-interactions (13) or, in the repul- 
sive case (18), at least to the scale |@| ~ Arep given in 
Eq. (23). The model should remain weakly coupled at 
these fields i.e. be renormalizable or have a sufficiently 
high cutoff. 


(b) The potential should be almost quadratic along 
the valley, V ~ m?|¢|?. All its repulsive terms should be 
suppressed at least by the scale Arep in Eqs. (18), (23). 
The potential may become flatter than quadratic (attrac- 
tive) inside finite ¢ intervals, but that should not ruin its 
renormalizability or destabilize the vacuum. If the po- 
tential becomes attractive asymptotically at || Z Aatt, 
like in Eq. (13), the scale Aatt should be sufficiently high, 
see Eq. (21). 

(c) Quantum corrections from the dark matter and vis- 
ible sectors should not destabilize the valley or create 
effective interactions breaking the condition (b). 

We have already demonstrated that the condition (c) 
is usually violated by dark matter interactions with the 
visible sector in one-field DM models. We now turn to 
models where this can be avoided. 


IV. MODELS WITH BENDED VALLEYS 
A. The mechanism 


Start from the arbitrary model for the DM field ¢ and 
add the second complex scalar x in such a way that ¢ 
and x have global charges 1 and 2, respectively. The 
coupling between the two is then chosen in a specific 
renormalizable and positive—definite form, 


V =d|d? — Ax] + m?o? +A |d|4/4, (30) 


where the last two terms represent the original @ po- 
tential. This model is invariant under phase rotations 
o — eh and y => e?"¥y. In the vacuum ¢ = y = 0 all 
fields are massive: mg =m and my = AA!/2. We will 
assume that the first term in the potential is the largest: 

m<AM? and N<dX. (31) 
This means, in particular, that the y—particles are heavy 
and do not change cosmology. 

The trick of the new model is to make the potential 
valley bend in the ¢—y space. Indeed, the low-energy 
field configurations are expected to minimize the largest 
term in Eq. (30), i.e. satisfy 


x= /A. (32) 


The other terms of V create a small potential along this 
valley. The same is true, in particular, for the stationary 
nonrelativistic soliton, 

p = p(x)",  x=ğ(æ)eTt, (33) 
which has real profiles y, ¥ satisfying Eq. (32). 

Let us explicitly demonstrate that the valley extends 
in the @ direction at small fields, then takes a turn at 
~ x~ A and goes along x. To this end we introduce 
a combination p(x) of the solitonic profiles which has a 


canonical kinetic term along the valley: (0;p)? = (3p)? + 
(0;X)?. Using Eq. (32) and integrating, we find 


p(y) = En X arcsinh (2) . (34) 


The classical potential V (p) at the bottom of the valley 
is obtained by inverting Eq. (34) and substituting y(p) 
into the last two terms of Eq. (30). 

At small fields we get p ~ y; hence, the valley stretches 
in the @ direction, indeed. Taylor series expansion in 
p/A <1 then gives a nonlinear valley potential, 


V(p) =m? p? + App*/4 + O/A) (35) 


with A, = X — 16m?/(3A?). We have already discussed 
in Sec. III that the pf term cannot be large positive, or 
it would stop the soliton from growing. Requiring 


N /4= B(m/A)’ , 


we ensure that the effective coupling is attractive at small 
fields: A, < 0. This is not dangerous for vacuum stabil- 
ity, since the overall potential of our model is explicitly 
positive—definite. 

The region ¢ > A is entirely different. Expression (34) 
gives p © y”/A meaning that the valley runs along ¥ ~% p. 
At the same time, the potential at the bottom of the 
valley is quadratic: V(p) = u?p?+(smaller terms), where 
we denoted 


Bots (36) 


p= mp? Sm, (37) 


and used Eqs. (30), (36). So, quartic self-interaction in 
this part of the valley is absent, and the effective p mass 
is smaller than m. 

Now we see how the bended valley (32) works. Orig- 
inally, the model (30) includes the term \|¢|* creating 
pressure. The same term, however, reduces to a mass 
|¢|* o |x|? once the valley takes a turn. Moreover, this 
property is valid even at the quantum level. 

Indeed, an explicit one-loop calculation [68] gives 
quartic effective potential at large fields: AVı—loop % 
Aa, eft |O|*/4 + O(¢?) with 


2 lel 
Ma, of = <z (20° + AV + 57/32) In 5 . (38) 


ren 


Thus, quantum fluctuations of ọ and x somewhat 
increase the constant A’ — A + A4, in front of 
Jolt œx |x|?, but do not prevent the Bose star from grow- 
ing if Eq. (36) remains satisfied: 


N + Aeg < 4(m/A)?. (39) 


The latter inequality is easily met if A is not too far away 
from m. 

One may wonder, why the dangerous terms AV = 
cılox|? + c2|x|4 are not generated at the one-loop level. 
If present, they would create an undesired pressure in- 
side the strong—field soliton. However, the auxiliary field 


x enters the interaction potential in the combination Ax, 
where A has mass dimension 1. On dimensional grounds, 
C2 ~ 2 m~ (A/Neutop)*, where Acutog is a cutoff for loops. 
Thus, the multiloop diagrams for the dangerous terms 
converge, the constants c1, C2 do not depend logarith- 
mically on the fields, and one can tune them to zero at 
o, xX > A. Barring this fine-tuning, the quantum correc- 
tions are harmless in the model (30). 

Let us visualize the growth of the Bose-Einstein con- 
densate in the model of this Section. Initially, it forms 
a Bose star held by the gravitational forces. This ob- 
ject becomes denser with time due to continuous inflow 
of dark matter particles. At p ~ m?/(|Ap|Mpi) < A the 
App“ self-attraction overcomes the kinetic pressure, and 
the Bose star collapses as a bosenova [54, 73]. This 
means that its central part starts squeezing in a partic- 
ular self-similar fashion [54, 55] developing strong fields. 
The squeeze—in halts when the field in the center reaches 
p~ A and the valley potential stops being attractive, 
cf. [55, 74]. 

This is the moment when a droplet of much denser 
condensate with p ~ x = A appears in the Bose star 
center. At first, it has the field strength y ~ A. How- 
ever, the droplet is expected to grow in density towards 
x >> A as more y-particles join in. Indeed, once it is 
there, a condensation process ¢¢ —> x with energy re- 
lease 2m — u becomes allowed. Even if one assumes that 
this direct process is ineffective dynamically, the conden- 
sation should continue in other, recurrent way. Without 
direct transmutation into x, the ¢ particles would cre- 
ate another dilute Bose star around the x—droplet, and 
that star would collapse, again, feeding the droplet. In 
any case all global charge should be eventually trans- 
ported into the dense y—-cloud and may never leave it 
due to large binding energy 2m — u of the x particles, cf. 
Eq. (37). 

Notably, we expect that the central y—cloud quickly 
thermalizes into a dense and non-interacting Bose star 
with x > A. Indeed, the y—particles effectively scat- 
ter with self-coupling = A’ at the cloud boundaries. In 
there, y ~ A and hence the particle number density is 
huge: n ~ pA?. A conservative estimate constrains the 
relaxation time from the above [43, 45, 75], 


2 100 GeV)? 
Ty ~ (oggun f) 1 K 40yr v —_ ; 


where € > (A/Mp)? is the time fraction spent by the 
particles at the cloud periphery, og, Z A? /(64ru?) is 
their cross section, f « n/(juv,)? is the phase-space den- 
sity, vy is velocity, and we used Eqs. (36), (37). Even 
the most conservative choices of parameters give small 
Ty compared to the lifetime of the Universe — thus, the 
relaxation is effective, indeed. 

Now, we recall that the self-interactions are absent at 
x >> A. This means that the respective Bose star will 
grow pressureless until it collapses gravitationally into a 
black hole. The latter event requires the critical charge 


Ner 21.3 Mj, /(8m") ’ (41) 


(40) 


where the numerical coefficient is larger by a factor 2/8 
than in Eq. (20) because x ~ p has global charge 2 and 
mass u = m61/?, 

We thus constructed a model with almost an optimal 
critical multiplicity for the collapse of Bose-Einstein con- 
densate into a black hole. In the next Section we will see 
that interaction with the visible sector and hence DM 
capture can be added to this model without disrupting 
the picture. 


B. Adding interactions with the visible sector 


We couple the DM field to the Higgs doublet in the 
same way as before — by adding the first term in Eq. (25) 
to the scalar potential (30). This makes the -particles 
scatter off neutrons with the cross section o x y? in 
Eq. (26). Importantly, we leave x to interact only with ¢. 
This maintains the desired properties of the strong—field 
Bose stars in our model. 

Now, the Higgs field changes along the potential val- 
ley (32) according to Eq. (27). As a consequence, a 
nonzero Higgs profile h(a) is generated inside every sta- 
tionary soliton. A combination of the three profiles — y, 
x, and h — with the canonical kinetic term is 


p Ay’? y? pl? 
=] dy’,/1 42 
p(y) f z t ua (42) 


where the last contribution comes from the Higgs field. 

Notably, with interactions (25) included, the effect of 
the bended valley remains essentially the same. Indeed, 
p is still proportional to y and ¥ at p< A and p > A, 
respectively. In the weak-field regime p < A the val- 
ley potential is nonlinear, Eq. (35), and has almost the 
same quartic constant A, as before. We therefore again 
impose the condition (36) to make Ap negative at weak 
fields. Once the valley turns at p > A, the potential 
V(p) becomes quadratic with mass (37), as guaranteed 
by the mechanism of the previous Section. 

All these nice properties remain valid because y does 
not couple directly to the Standard Model fields and still 
enters the overall potential in the combination Ay with 
dimensionful A. This makes the nonlinear interactions 
vanish in the strong—field region, and they cannot be gen- 
erated by loops. We check the latter property by com- 
puting the one-loop effective potential in the three—field 
model of ¢, x, and H at the bottom of the valley (27), 
(32). The result at large fields is given by Eq. (28), where 
the constant 4, e equals Eq. (38) plus the Higgs contri- 
bution 


2 j : 
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Similarly, the other Standard Model fields are expected 
to generate AVi—loop « h* œ ||*. Regardless of the 
value of A4 ef, the loop corrections safely change the mass 
||* x |x|? in the strong-field region. The value of u then 


can be made positive by adjusting A’. For simplicity we 
assume below that the loop corrections are small com- 
pared to A. 

Without the interaction pressure, the strong—field Bose 
star collapses gravitationally at almost optimal critical 
charge (41). This value cannot exceed the total number 
of captured dark matter particles, More /m > Ner. 
One obtains the inequality 


m 1/2 a 
100 ay) B =) 
using Eqs. (6), (26), and (41). This constraint is rela- 
tively mild and can be easily satisfied. 

Indeed, let us show that a viable parametric region for 
our model can be selected. One starts by specifying m 
within the WIMP mass range and 8 < 1. The coupling 
constant y is then chosen to satisfy Eq. (43) together 
with the constraints from the DM detection experiments. 
Below we will argue that a narrower region 


y > 1077 ( 


m 


hp <2x 10-5 671/2 (44) 
100 GeV È 


2x 10-8 sal 


is better for phenomenology. The next parameters are A 
and \’ < à. In the simplest case one takes À ~ y/ A 
and X’ = A?, thus suppressing the running of \’: A4, ef, 
d4, et SX’. Finally, the scale A is given by Eq. (36). 
We conclude that all the constraints can be satisfied if 
the hierarchy 

XN ~X wy Àn, 


Anm/X, BSl1 (45) 


is valid and y belongs to the interval (44). 


C. Generalizations 


So far we relied on fine-tuning to kill the |¢x|? and 
|x|* terms in the potential. If added, these interactions 
should be extremely weak — say, the operator \4|y|*/4 
should satisfy the analog of Eq. (24): 


1 < 2x 107! $263 (m/100 GeV)* . 


Can we allow stronger interactions at |x| Z A? 

The answer is yes, if extra auxiliary scalars are intro- 
duced. Namely, let us add the |x|* term along with the 
new field x2 that has a global charge 4. The new terms 
in the potential are 


Va = Aalx? — Aaxal? + ASlxI*/4, (46) 


where A5 = 462(u/A2)? and b2 < 1; cf. Eq. (36). Now, 
the valley takes a second turn at x ~ x2 ~ Ag, before the 
||* interaction becomes relevant. 

One can continue this procedure and arrive at the 
“clockwork—like” model with n fields y;, a hierar- 


chy of scales A < Ag < < An, masses at 
strong fields Hi = Hi—1 8, / . and many coupling constants 


A; = (2u;/A;)?, where 8; < 1. The Bose star made of 
the last field Xn collapses gravitationally at the critical 
charge 
n 
Nor = (Myi/m)? |] (2/6:) - (47) 
i=1 
If n ~ O(1), this multiplicity is not exceedingly large. 
Of course, the field y can be made interacting in other, 
old-fashioned ways — say, by supersymmetry. To this 
end one upgrades this field to the supersymmetric sector 
with a flat direction [76] which represents the |x| Z A 
part of the valley. The flat direction should be lifted by 
soft terms giving the mass to the condensate. Then the 
superpotential will be protected from quantum correc- 
tions as long as its fields couple to the visible matter via 
the dimensionful constant A. 


V. CONCLUSIONS 


We considered formation of black holes with 
masses © 1 Mo by means of bosonic dark matter collapse 
inside neutron stars. This scenario includes DM capture 
by the neutron stars, its thermalization with neutrons, 
Bose-Einstein condensation, and at last — gravitational 
collapse into seed black holes which eventually consume 
the stars. The overall process includes two conflicting re- 
quirements on the DM model [18]. On the one hand, the 
DM can be captured only if it interacts strongly enough 
with the visible sector. On the other hand, loop con- 
tributions of the same interactions generate DM self- 
couplings and hence pressure inside dense DM clouds, 
impeding their collapse. We have found, in agreement 
with Ref. [18], that the conflict cannot be resolved by 
tuning the parameters of the single-field DM models, 
neither can it be achieved by optimization of their scalar 
potentials. 

Using the crucial observation that the conflicting re- 
quirements are imposed at different scales, we proposed 
a mechanism that enables black hole formation within 
this scenario. The respective DM models are deformed 
at strong fields, i.e. in the regions inaccessible to the di- 
rect experiments. Their scalar potentials include bended 
valleys that go along the dark matter scalar ¢ at |¢| S A 
and then turn in the direction of a new field y. In this 
case x self-pressure can be made small and, notably, 
can be protected from all quantum corrections by super- 
renormalizable couplings to other fields, supersymmetry, 
or a “clockwork-—like” mechanism. As a consequence, the 
growth of the DM cloud inside the neutron star includes 
a phase transition ġġ —> x at a certain density followed 
by the gravitational collapse of the pressureless x con- 
densate. 

In our model formation of a black hole requires a rela- 
tively small number Eq. (41) of DM particles which can 
be accumulated inside the neutron star even at extremely 
weak couplings, see e.g. Eq. (43) in the case of interac- 
tion (25). However, other conditions become relevant at 
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this point. First, the newborn seed black hole cannot 
be too light, Mgy = 101% g, or it evaporates faster than 
it accretes neutrons [16]. In our model Mpy ~% uNer/2, 
so this condition constrains an effective mass u of the 
auxiliary scalar x at strong fields, 


u = mB? < 17 GeV, (48) 


where Eq. (41) was used, m is the mass of the DM par- 
ticles, and 8 <1 parametrizes their self—interactions via 
Eq. (36). Taking sufficiently small 6, one fulfills Eq. (48) 
in the entire WIMP mass range. 

Second, thermalization of DM with neutrons is an 
obligatory part of our scenario needed to form a dense 
central cloud. But at weak couplings the equilibrium 
may be unachievable even on the cosmological timescales, 
since DM interactions with neutrons are further sup- 
pressed at low energies by Pauli blocking [14]. For scalar 
vertices, the respective thermalization time was esti- 
mated in [42] as tus © 357? /(12omT?), where T ~ 10° K 
is the neutron star temperature and ø is the DM- 
neutron cross section in vacuum, Eq. (26). Requiring 
Tys Í 101° yr, we constrain the coupling constant, 

1/2 
m ) / (49) 


>72 10 (y 
MA | y 


This condition is marginally stronger than Eq. (43) 
needed for capture. 

Third and finally, our mechanism should not be too ef- 
ficient in turning all neutron stars into solar—mass black 
holes — after all, thousands of neutron stars are observed 
in our Galaxy. One way to suppress their transmutation 
is to make the DM-neutron coupling moderately small, 
so that only the densest DM environments would allow 
the stars to accumulate the critical mass for collapse. Us- 
ing the Milky Way parameters, we require MMW < MNer 
or 


m 12 =1/2 
100 = P , (50) 
where Eqs. (5), (26), and (41) were used. Note that 
Eq. (50) is not far from the conditions (43), (49); this 
is due to the fact that the DM parameters vary by only 
a few orders of magnitude from galaxy to galaxy. 

It is worth noting that we have disregarded several 
points which are important for the scenario of this paper. 
First, it is crucial that the DM is non-annihilating and 
satisfies constraints coming from its generation in the 
early Universe and from the DM detection experiments. 
It remains to be seen if all these conditions can be fitted 
together with the requirements of our mechanism. 

Second, we have not analyzed the astrophysical signa- 
tures of the neutron stars converting into the solar—mass 
black holes. In our scenario, the transmutation is con- 
trolled by a single parameter — the number of DM par- 
ticles Ner required for collapse. Thus, regardless of the 
underlying DM model the neutron stars will be converted 
in the parts of the Universe where the DM abundance 


yS2x 10-5 ( 


and velocity ensure its efficient accumulation. As a con- 
sequence, solar—mass black holes will be distributed in a 
specific way among the galaxies of different types. For 
example, a discovery of old neutron stars in the dwarf 
galaxies — the best-known environments for the DM ac- 
cumulation — would exclude a sizable overall abundance 
of the transmuted neutron stars in the Universe. 

Third, several parts of our mechanism deserve a de- 
tailed numerical study. We expect that the phase tran- 
sition of the DM particles into the x quanta proceeds in 
a spectacular first-order way, starting as a self-similar 
bosenova collapse [54, 55] and ending up with formation 
of a dense y condensate. This full two-stage proscess 
has never been simulated before; in the main text we 
have just crudely estimated the time of x condensation 
to be small, see Eq. (40). The other unexplored subject 
is the growth of the y Bose-Einstein condensate in the 
end of the process. We assumed that the growth contin- 
ues even when the condensate becomes incredibly small 
in size. On the one hand, this optimism is based on the 
fact that unlike in Refs. [45, 46, 77, 78], the DM interac- 
tions in our model are short-range. On the other hand, 
even in the worst case the growth of the -condensate 
may proceed in a recurring indirect way: by growing the 
@ condensate via thermalization to the point when the 
transition ¢¢ > x happens, and then growing the ¢ con- 
densate, again. The details of this process also should be 
studied numerically. 

Fourth, although viability of our new mechanism does 
not depend on the details of dark matter capture and 
thermalization, its predictions in a particular model rely 
on a quantitative description of these phenomena. Thus, 
precise identification of a phenomenologically acceptable 
parameter region for the dark matter model should in- 
clude studies of astrophysical uncertainties, equation of 
state for the nuclear matter, and DM interactions with 
it. Besides, renormalization group equations should be 
solved for the dark matter constants \’ and A which are 
scale—dependent and run on par with the Standard Model 
couplings. 

Fifth and finally, it might be interesting to explore 
dark matter collapse in the centers of white dwarfs within 
our model. Despite being more dilute than the neutron 
stars, the latter objects are also better understood, which 
makes them a perspective testing ground for dark matter 
studies. 

In a nutshell, our study demonstrates that almost any 
model of bosonic DM can be modified at strong fields in 
such a way that the solar-mass black holes can appear 
by transmuting the neutron stars. This calls for a proper 
identification of low-mass compact objects in astrophys- 
ical observations, cf. [79] and [80]. 
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Appendix A: Q—ball and Bose star parameters 


Here we derive parametric dependence of the soliton 
mass M and radius R on its global charge N in the mod- 
els of Sec. III. To this end we adopt a simple version of 
the variational ansatz [58, 73, 81] ignoring order—one nu- 
merical coefficients wherever possible. We will separately 
study the solitonic Q—ball [82] stabilized by the attrac- 
tive self-interactions and the gravitationally bound Bose 
star [83] in the case of self-repulsion. 

Start with the Q-ball in the model (9) for the com- 
plex scalar field ¢(a). We will assume that its scalar 
potential V becomes flat, i.e. attractive, at strong fields 
|d| Z Aatt, where it can be roughly approximated as a 
power-law (13) with a < 2. Suppose the stationary Q- 
ball in this model has the form of a single bell-shaped 
lump with typical field strength yo and size R: 


$ = of(x/R)e™ , 


where f(y) is a dimensionless order—one function with 
the support at |y| = |a|/R < 1; the time dependence fol- 
lows from Eq. (11). 

Since the Q-ball is mostly bound by the self- 
attraction, we can compute its parameters in flat space- 
time ignoring gravity. Substituting Eq. (A1) into the 
energy and charge of this object, we get, 


(Al) 


Mw wpe R? + yR + mye ih2 AR i 
Nw wpe R? ; 


(A2) 
(A3) 


where we omitted yo- and R-independent numerical co- 
efficients of order 1 in front of every term — these come 
from the integrals over y of powers of f(y) and its deriva- 
tives. 

We now minimize M with respect to yo and R at a 
fixed N. Expressing w from Eq. (A3), we substitute it 
into Eq. (A2) and differentiate the result with respect to 
the unknowns. We get, 


ƏM 2N? aaa 
£0 Ovo ms y2 R3 2p R t am po Au R? =0 ’ 
0 
OM 3N2 ae 
R IR ~~ ake PÈR + 3m pth? R? =0. 
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Notably, the solution of these equations does not exist 
at a > 2 when self—attraction changes to self-repulsion. 


Away from that region, all terms in the equations are of 
the same order: 


oo ~ (mR)? AGS and N~ yo R? ; (A4) 


Using Eq. (A3), we find that w ~ R7! and reproduce 
Eqs. (14), (15) from the main text. 


It is worth noting that w ~ R7! < m at sufficiently 
large N, see Eq. (15). Thus, the particles inside the large 
Q-ball are very light due to mass deficit. As a conse- 
quence, the overall soliton is lighter and less amenable to 
gravitational collapse than a collection of free particles 
with the same charge. 


Now, consider a Bose star in the model with the scalar 
potential (18) and a > 2. Recall that the flat—space solu- 
tion does not exist in this case — hence, we add gravita- 
tional attraction. In the stationary non-relativistic limit, 
this amounts to introducing the metric goo = 1+2U and 
gij = —(1 — 2U) ĝij, where U < 0 is a time-independent 
Newtonian potential, |U| <1. With gravity included, 
the energy and charge of the nonrelativistic soliton read, 


M= f a {(1—0)lasl? + [a6 + (U)?/(8x6) 
+ (1= 2U) m?l? + mèla}, (A5) 


wai J dæ (1— 4U)($*dod — pog"), (A6) 


where the last term in the first line represents the energy 
density of the gravitational field and we ignored the grav- 
itational effect of the self-interaction energy ôV « |¢|®. 


Following the same strategy as before, we assume that 
the solution has the form (A1) characterized by the field 
strength yo, size R, and the gravitational potential in the 
center Up. Omitting all order—one coefficients, again, we 
find, 


M ~ (1—4Uo) w? pi R? + pR +UR/G 
+ (1 — 2U0) mpi R? + i el R? ; 
N ~ (1—4Up) wy R? , 


cf. Eqs. (A2), (A3). Below we also disregard the gradient 
energy of the soliton y2R, checking a posteriori that it is 
small. 


Next, we express w from Eq. (A8), substitute it into 
Eq. (A7), and minimize the resulting energy with respect 
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to yo, R, and Uo. The extremality conditions are, 


(1 + 4U0) ae + 2(1 — 2U0) m pR? 
: + am ph AZo RS =0, (A9) 

(1+ 40) oF + 3(1 — 2Uo) m? ye R? 
+ US R/G + 3mp At R? =0, (A10) 
D + aon -mp R? =0, (A11) 


where we recalled that |Uo| « 1. It is peculiar that the 
above system is almost degenerate: Eqs. (A9) and (A10) 
differ only by the last terms suppressed as Uj or A2. 
This is a consequence of the fact that the dominant part 
of the potential is quadratic. To the leading order, both 
equations give the standard expression for the number of 
nonrelativistic particles, 


N ~ mR? . (A12) 

Now, Eq. (A11) takes a familiar form Up ~ —GmN/R. 

The last relation is found from the linear combination of 

Eqs. (A9) and (A10) such that the large terms cancel. 
We get, 

2UG/G ~ 3 (a — 2) (mR) pe A27*, 


rep 


(A13) 


Since the left-hand side is positive, the solution exists 
precisely in the case aœ > 2. One finally obtains, 


a—4 
Ẹ 3a—8 j 


2 Mpy 
~ A. 3a—8 Rw pair a 
Yo rep E ry Arep 


(A14) 


where € = Nm? Arep/M3,. 

The interpretation of the above solution is essentially 
different in the cases 2 < a < 8/3 and a> 8/3. Ifa 
exceeds 8/3, the Bose star is stable at € < 1. Its field 
steadily grows with the multiplicity reaching yo ~ Arep 
at E ~ 1. At this point, the soliton’s gravitational po- 
tential Up ~ —€2—-4)/GBe-8) becomes of order one and 
a black hole forms, see Eq. (A13). Notably, at the brink 
of collapse mR ~ Myi/Arep > 1 — thus, the soliton 
remains nonrelativistic, indeed. Also, the gradient term 
pR in Eq. (A7) is much smaller than the binding energy 
U§R/G, as we assumed. 

In the opposite case 2 < a < 8/3 the scalar self- 
repulsion triumphs over gravity at large distances. In- 
deed, expressing yo from Eq. (A12), one learns that at 
a fixed N the self-repulsion energy of the nonrelativistic 
condensate grows faster with its size R than the grav- 
itational binding energy: (y¢R°)/(UZR) œx R873%)/2, 
see Eq. (A7). This means that clumps of Bose-Einstein 
condensate spread over the entire volume offered by the 
external conditions. The solutions (A14) in this case 
break the Vakhitov—-Kolokolov criterion [54, 84] dw/dN ~ 
mdUo/dN < 0 necessary for stability. Thus, they repre- 
sent the maxima of the potential energy and determine 
the minimum size R to which the condensate should be 
squeezed for gravity to dominate. We do not consider 
any external forces in this paper — hence, the black hole 
does not form in the case 2 < a < 8/3 at all. 
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